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ABSTRACT 

Oil  consumption  was  measured  in  a  modem  production  diesel  engine  using  tritium  as  a 
radiotracer.  The  measurements  were  made  primarily  at  two  speeds  and  one  load  using  first 
a  single-grade  lubricant  and  then  a  multi-grade  lubricant. 

These  values  were  then  compared  to  oil  flow  rates  up/down  the  liner  which  were  based 
on  film  thickness  traces  of  a  sister  engine  under  the  same  loads  and  speeds.  The  traces 
were  obtained  using  the  laser-fluorescence  technique.  For  the  most  part,  it  was  discovered 
that  there  does  not  seem  to  exist  a  correlation  between  these  flow  rates  and  oil 
consumption.  However,  the  traces  do  reveal  that  the  crown  land  is  dry  on  all  four  strokes 
and  thus  does  not  contribute  to  the  engine’s  oil  consumption. 

A  larger  data  base  is  necessary  in  order  to  :curately  compare  oil  consumption  to  the 
film  traces.  This  is  currently  in  progress  as  of  this  writing.  f 
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CHAPTER  1  BACKGROUND  AND  INTRODUCTION 

Due  to  the  EPA  particulate  emission  standards  for  1991  and  1994,  it  has  become 
necessary  to  learn  more  about  oil  consumption  and  the  driving  forces  behind  it. 
Subsequently,  an  accurate  determination  of  these  driving  forces  could  lead  to  designs  of 
pistons,  rings,  and  lubricants  which  improve  both  oil  consumption  and  paniculate 
emission. 

This  paper,  then,  has  two  purposes.  First,  initiate  an  industry  proven  method  for 
measuring  oil  consumption,  and  second,  compare  these  values  to  the  average  oil  flow  rates 
up  the  liner  of  a  sister  engine  operating  at  the  same  loads  and  speeds.  The  oil  consumption 
measurement  technique,  along  with  engine  >et-up  and  experimental  procedures,  will  be 
explained  further.  The  oil  flow  rates,  however,  have  been  calculated  from  oil  film 
thickness  data  taken  from,  and  explained  in  separate  papers  [1,2  respectively]. 
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CHAPTER  2  DESCRIPTION  .OF  MEASURING  TECHNIQUE 

The  technique  used  is  a  radiotracer  method  which  uses  tritium  as  its  radiotracer.  This 
method  has  been  widely  used  for  assessing  oil  consumption  for  over  two  decades.  The 
method  is  fast,  accurate,  and  economical  relative  to  non-radioactive  methods.  This  overall 
practicality  and  effectiveness  is  a  major  result  of  the  advances  of  liquid  scintillation 
counting. 

2.1  BASIC  METHOD 

The  following  is  a  brief  description  of  the  basic  method  while  a  more  detailed 
description  of  the  general  technique  is  found  in  Reference  [3],  However,  for  ease  of 
reference,  the  derivation  of  the  oil  consumpion  calculation  equation,  upon  which  this 
method  is  based,  is  included  as  Appendix  A. 

The  basic  procedure  consists  of  running  an  engine  with  radioactive  lubricating  oil  with 
a  known  radioactive  level  of  disintegrations  per  minute  (dpm).  The  total  amount  of 
hydrogen  passing  through  the  engine,  via  fuel  flow  and  air  flow,  is  measured  and 
calculated  Through  combustion,  this  hydrogen,  along  with  the  hydrogen/tritium  from  the 
oil.  leads  to  the  formation  of  water.  A  sample  of  this  water  is  collected  by  condensing  a 
sample  of  the  exhaust  gas  (Only  a  sample  is  needed  since  total  measurements  of  the  flow 
rates  into  the  engine  are  obtained)  Subsequently,  the  radiation  concentrations  of  both  the 
condensed  water  and  the  oil  are  detennined  using  a  liquid  scintillation  counter.  These 
values  are  then  applied  to  the  oil  consumption  equation  to  arrive  at  the  engine's  total  oil 
consumption 

2.2  SOOT  SAMPLE  COLLECTION 

As  seen  in  Appendix  A.  the  general  technique  is  based  on  complete  combustion  of  the 
exhaust  sample.  For  this  paper,  however,  the  exhaust  sample  is  Filtered  and  then 
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condensed  instead  of  being  completely  combusted  and  then  condensed.  This  results  in  an 
activity  concentration  (dpm/ml)  of  the  condensed  water  and  a  total  activity  (dpm)  of  the  soot 
sample. 

In  order  for  this  method  to  be  valid,  two  assumptions  are  made.  First,  the  soot 
collected,  if  burned,  would  produce  a  negligible  amount  of  water.  Secondly,  all  of  the 
water  vapor  in  the  sample  is  being  condensed.  This  latter  assumption  is  necessary  since  the 
consumption  calculation  requires  radioactive  concentrations  vice  total  activities.  Therefore, 
the  soot  sample  total  activity  is  divided  by  the  amount  of  water  collected  and  then  added  to 
the  activity  concentration  of  the  water  sample.  A  more  in  depth  description  of  the 
procedures  is  included  in  Appendix  E. 

2.3  LIQUID  SCINTILLATION  COUNTING  f41 

All  radioactive  levels  are  determined  by  a  method  known  as  liquid  scintillation  counting. 
This  technique  is  highly  sensitive  and  greatly  responsible  for  the  increased  use  of  weak  beta 
emitters  such  as  carbon- 14  and  tritium.  The  basic  technique,  as  with  other  methods  of 
measuring  radioactivity,  relies  on  the  interaction  of  nuclear  radiation  with  matter.  In  this 
case,  a  luminescent  material,  or  scintillator,  is  excited  by  radioactive  disintegrations. 

A  scintillator,  also  referred  to  as  either  a  phosphor  or  sometimes  a  fluor,  is  defined  as  a 
materia]  which  emits  a  brief  pulse  of  fluorescent  light  when  interacting  with  a  high-energy 
particle  or  quantum.  Therefore,  the  radioactive  samples  used  in  this  experiment  are  mixed 
with  a  scintillator.  Tire  subsequent  flashes  of  light,  or  scintillations,  are  detected  by  the 
counter  by  use  of  a  photomultiplier.  This  in  turn,  over  a  set  period  of  time,  allows  for  the 
determination  of  the  average  disintegrations  per  minute  (dpm)  of  each  sample. 


CHAPTER  3  E 


Ol'IPMENT  AND  SET-LIP 


3.1  ENGINE 

The  engine  used  in  this  experiment  is  the  KUBOTA  EA300N,  a  single  cylinder, 4- 
stroke.  ID1  diesel.  This  engine  is  used  primarily  for  remote  power  generation.  A  detailed 
description  of  the  engine's  geometry  and  performance  is  included  in  Appendix  B. 

This  engine  was  chosen  only  because  the  film  thickness  measurements  used  in  this 
paper  were  obtained  from  an  engine  of  the  same  make  and  model  [1  ]. 

3.2  RADIOACTIVE  LUBRICATING  OILS 

Radioactive  oil  for  these  experiments  is  purchased  commercially.  The  oil  received  has 
been  through  a  tritiation  process  which  results  in  approximately  250  mCu  in  one  ml  of  oil. 
Subsequently,  this  oil  is  divided  up  and  thoroughly  mixed  with  Pennzoil  SAE-30  (single¬ 
grade)  and  Pennzoil  15W-40  (multi-grade).  The  tritiation  process  is  explained  further  in 
Appendix  B. 

These  oils  are  used  primarily  for  the  same  reason  as  above  for  the  engine.  In  addition, 
they  should  provide  a  good  comparison  of  single-  and  multi-grade  oils  since  the  only 
difference  is  the  addition  of  a  polymer  VI  improver  to  the  1 5W-40.[3] 

3.3  ENGINE  INSTRUMENTATION 

The  engine  is  fitted  with  thennocouples  n  order  to  monitor  operating  temperatures  (Fig. 
1 ).  The  thennocouples  measure  the  inlet  air  temperature,  the  coolant  temperature,  the  oil 
sump  temperature,  the  exhaust  temperature,  and  the  liner  temperatures.  The  latter 
thermocouples  coincide  with  the  top  ring  at  TDC.  mid-stroke,  and  BDC.  A  laminar  air 
flow  element  is  also  used  to  measure  the  air  flow’  while  a  burette  system  is  used  to  measure 
the  volume  flow  rate  of  the  fuel  A  more  detailed  listing  of  the  measuring  equipment  is 
included  in  Appendix  B. 
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3.4  DESCRIPTION  OF  PISTON  AND  RINGS  121 

The  piston  and  rings  (Fig.  2)  are  representative  of  modem  engine  design.  The  piston  is 
made  of  aluminum  and  accommodate  three  piston  rings.  It  has  round  lands  and  a  barrel 
shaped  skirt.  Behind  and  below  the  oil  control  ring,  there  are  oil  flow  relief  holes  that 
allow  the  circulation  of  oil  up  the  piston  skirt  on  the  liner  side  and  back  down  to  the  sump 
inside  the  piston. 

The  three  piston  rings  consist  of  a  top  ring,  a  scraper  ring,  and  an  oil  control  ring.  The 
top  ring  is  a  half-keystone  design  with  a  7°  keystone  angle.  It  is  made  out  of  cast  iron  with 
a  chromium  strip  embedded  into  its  face  to  improve  wear  resistance.  The  resulting  face 
profile  is  asymmetric.  The  scraper  ring  has  a  tapered  face  and  an  undercut  on  its  face  side 
to  provide  the  correct  tilt  angle  when  compiessed  to  its  working  diameter.  It  is  also  made 
of  cast  iron.  The  oil  control  ring  is  a  one  pi?ce  type  made  of  cast  iron  also,  and,  like  the  top 
ring,  its  rails  have  a  chromium  strip  embedded  into  them.  The  oil  control  ring  also  has 
slots  between  the  rails  to  allow  the  flow  of  oil  through  the  ring  to  the  holes  in  the  piston 
mentioned  above. 

3.5  SAMPLING  SYSTEM 

Figure  3  is  a  schematic  of  the  sampling  system  in  use.  As  shown,  the  sampling  system 
is  connected  to  an  exhaust  tank.  This  tank  is  used  to  help  eliminate  the  exhaust  pulse  that 
exists  when  operating  a  single  cylinder  eng'ne.  This  entire  system,  consisting  of 
inexpensive  components,  is  placed  on  a  small  lab  cart  which  allows  for  portability  and  for 
an  easy  connection  to  a  different  engine. 
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3.6  OIL  FILM  THICKNESS  MEASUREMENTS 

The  film  thickness  traces  referred  to  are  obtained  from  reference  [1]  which  uses  a  laser- 
fluorescence  technique.  This  basic  technique  consists  of  mixing  the  lubricant  with  a 
fluorescent  dye  and  then,  during  operation,  focusing  a  HeCd  (blue)  laser  through  a  quartz 
window  installed  in  the  liner.  The  resulting  fluoresced  (green)  light  is  then  collected  and 
converted  to  a  voltage  signal  using  a  photomultiplier  [2],  Lastly,  the  voltage  signal  is 
converted  into  microns  using  a  calibration  coeffecient  (unique  for  each  trace).  This 
coeffecient  is  based  on  the  voltage  readings  obtained  when  etch  marks  of  known  depths  on 
the  piston  skirt  pass  the  quartz  window.[l] 

More  in  depth  descriptions  of  this  technique  can  be  found  in  references  [1]  and  [2]. 
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CHAPTER  4  DESCRIPTION  OF  EXPERIMENTS 

4.1  PROCEDI  RES 

At  each  given  load  and  speed,  the  engine  is  allowed  to  stabilize.  Subsequently,  flow 
rates  are  measured  and  samples  are  taken  in  accordance  with  the  step-by-step  procedures  in 
Appendix  C.  The  samples  are  then  prepared  in  accordance  with  Appendix  D  in  order  for 
the  scintillation  counter  to  be  effective. 

Following  sample  preparation,  the  samples  are  placed  in  the  scintillation  counter  to 
determine  their  individual  radioactive  levels.  These  values,  along  with  the  measured  flow 
rates,  are  then  used  to  calculate  oil  consumption.  Actual  formulas  and  an  example  of  the 
spreadsheets  in  use  are  included  in  Appendix  E. 

After  collecting  data  using  the  first  oil,  the  engine  is  flushed  to  minimize  contamination 
of  the  next  test  oil.  A  flow  chart  of  these  basic  procedures  is  illustrated  in  figure  4. 

4.2  ENGINE  OPERATING  CONDITIONS 

The  primary  loading  conditions  for  each  oil  consisted  of  full  load  at  1 500  rpm  and 
full  load  at  3000  rpm  (approximately  9.5  ft-lbs  for  each).  As  before,  these  conditions  were 
used  in  order  to  duplicate  the  conditions  under  which  the  film  thickness  data  was  taken. 
Secondary  loading  conditions  consisted  of  2000  rpm  and  2500  rpm  at  the  same  torque. 
These  latter  conditions  were  used  mainly  fc  r  graphic  continuity  between  low  and  high 
speeds. 
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CHAPTER  5  GRAPHIC  RESULTS 

5.1  OIL  CONSUMPTION  MEASUREMENTS 

Oil  consumption  measurements  for  the  desired  load  and  speeds  using  Pennzoil's  SAE- 
30  and  15W-40  are  illustrated  in  figures  5  and  6,  respectively.  Figure  7  is  then  an  overlay 
of  the  previous  two  graphs  with  figure  8  being  another  comparison  but  in  the  form  of  a  bar 
chart. 

The  volatility  of  a  new  oil  (SAE-30)  is  also  examined.  The  results,  included  as  figure  9, 
illustrate  how  oil  consumption  rates  decrease  and  then  stabilize  with  time  (for  a  constant 
load  and  speed).  Volatility  information  is  desirable  since  the  film  thicknesses  were 
obtained  from  the  sister  Kubota  using  new  lubricating  oils. 

5.2  OIL  FILM  THICKNESSES 

Figures  10.  11,  12,  and  13  illustrate  a  subset  of  the  average  film  thickness  traces  for  ten 
compression,  expansion,  exhaust,  and  intake  strokes,  respectively  [1],  For  these  traces, 
the  engine  was  using  SAE-30  and  was  running  at  full  load  at  1500  rpm.  The  traces  for 
3000  rpm  and  for  the  different  oil  are  similar  to  these  and  are  therefore  not  included. 
Interpretation  of  this  data,  which  can  be  confusing  to  the  first-time  observer,  is  described  in 
length  in  reference  [2]. 


S 
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CHAPTER  6  ANALYSIS  AM)  DISCUSSION 

6.1  HYPOTHESES 

There  exists  three  prevalent  hypotheses  on  lubricating  oil  consumption,  two  of  which 
deal  with  oil  being  drawn  in  and  around  the  top  ring  at  or  near  the  end  of  the  expansion 
stroke.  At  this  point,  the  high  pressure  below  the  top  ring  forces  the  oil  around  the  top  ring 
where  it  either  (1)  accumulates  on  the  crown  land  from  where  it  is  later  burned  off  or  (2)  is 
"blown”  up  into  the  combustion  chamber  and  consumed.  The  third  hypothesis  is  based  on 
oil  being  transported  by  the  top  ring  face.  Subsequently,  some  of  this  oil  is  deposited  onto 
the  liner  at  or  near  top  ring  reversal  (TRR)  where  it  is  later  consumed. 

For  each  hypothesis,  there  exists  a  corresponding  volume.  These  volumes  are 
calculated  by  first  finding  their  respective  aieas  which  are  illustrated  in  figure  14.  The  areas 
are  then  swept  around  the  circumference  of  the  bore  to  determine  the  volumes.  Table  1  lists 
the  areas  and  the  volumes  calculated. 


SPACE 

AREA  (mm2) 

VOLUME  (mm-) 

A)  Crown  land  crevice 

1.33 

313.0 

B)  Ring/groove  crevice 

0.3 

70.7 

(top  and  bottom) 

C)  Ring/groove  crevice 

1.02 

241.0 

(backside  of  ring) 

D)  Between  ring  and  liner 

.00386 

0.91 

Table  1.  Areas/Volumes  of  possible  oil  consumption  mechanisms 
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These  volumes  can  in  turn  be  related  to  an  oil  flow  rate,  in  g/hr,  based  on  engine 
operating  speed  and  oil  density  (Refer  back  to  Appendix  E).  Unfortunately,  volumes  A,  B, 
and  C  above  result  in  rates  on  the  order  of  thousands  of  grams  per  hour  with  volume  D 
being  on  the  order  of  50  grams  per  hour.  These  values  are  extraordinarily  high  and  do 
not  scale  with  oil  consumption  (figures  5  ard  6).  Thus  they  cannot  offer  much  guidance  in 
determining  oil  consumption  mechanisms.  Actual  volumes  necessary  to  resemble  the  oil 
consumption  values  are  on  the  order  of  0.0*%  of  volumes  A  and  C,  0.05%  of  volume  B, 
and  just  3.0%  of  volume  D. 

6.2  ANALYSIS  AND  DISCUSSION 

As  discussed  in  Chapter  5  of  reference  [2],  the  film  thickness  data  does  not  usually 
allow  for  a  determination  of  where  the  oil  lies  (i.e.  on  the  piston,  on  the  liner,  or  on  both). 
However,  as  figures  10  through  13  illustrate,  the  oil  film  thickness  before  the  piston  arrives 
(-4mm  to  Omm  for  the  upstrokes)  is  relatively  the  same  as  when  the  crown  land  is  passing 
(0mm  to  9.5mm  for  the  upstrokes).  Therefore,  it  is  apparent  that  the  piston  crown  land  is 
virtually  oil  free  and  thus,  does  not  contribute  to  oil  consumption.  It  is  important  to 
mention  that  this  film  trace  characteristic  is  also  shared  by  the  other  traces  for  the  higher 
speed  and  with  the  different  oil. 

Based  on  the  above  line  of  reasoning,  the  first  hypothesis  can  be  easily  refuted.  The 
second  one,  however,  can  neither  be  directly  refuted  nor  explicitly  proved  by  this  method 
of  film  measurement.  This  is  because  the  hypothesis  addresses  oil  which  enters  the 
combustion  chamber  by  way  of  neither  the  piston  nor  the  liner.  However,  further  data 
analysis  of  the  film  traces  does  allow  for  some  broad  conclusions  to  be  made  concerning 
both  this  and  the  final  hypothesis. 

This  subsequent  data  analysis  first  requires  the  use  of  existing  computer  programs 
which  have  been  modified  for  the  Kubota  engine  [5].  These  programs  calculate  the  volume 
differences  (mm-Vcycle)  of  various  crank  angle  regions  for  both  the  power  and  the  gas 
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exchange  strokes.  Only  five  regions,  which  are  listed  below  in  Table  2,  are  analyzed  in 
this  paper.  The  methods  and  reasons  for  computing  these  particular  regions,  along  with 
copies  of  the  programs,  are  included  in  Appendix  F. 


REGION 

PISTON  AREA 

I 

Crown  land 

11 

Top  Ring 

in 

1  ring  width  below  top  ring 

IV 

Second  Land 

V 

Second  (Scraper)  Ring 

Table  2.  Piston  Regions  of  Film  Thickness 


The  title  of  each  region  in  the  table  explains  which  part  of  the  piston  is  passing  the 
quartz  window  at  that  instant.  In  addition,  the  regions  are  labeled  on  an  expanded  version 
of  the  aforementioned  compression  stroke.  (Figures  15  and  16) 

Secondly,  in  order  to  compare  a  regions  volume  difference  to  the  measured  oil 
consumption,  the  former  is  converted  to  an  oil  flow  rate  up  (or  down)  the  liner  in  g/hr. 
This  is  done  based  on  the  volume  difference,  the  speed,  and  the  oil's  density.  (Refer  back 
to  Appendix  E  for  calculations) 

Tables  3  and  4  and  Tables  5  and  6  tabulate  the  data  collected  using  SAE-30  and  15W- 
40,  respectively.  A  positive  value  indicates  that  the  upstroke  volume  is  greater  than  the 
downstroke  and  vice  versa.  Subsequently,  figures  17  through  22  are  included  to  illustrate 
these  various  flow'  rates  of  the  two  oils  for  the  different  speeds.  On  each  figure,  the 
respective  oil  consumption  is  included  in  order  to  make  relative  comparisons. 
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PWR  EXCHANGE 

GAS  EXCHANGE 

TOTAL  (pwr+gas) 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

I 

0.04 

1.51 

0.00 

0.00 

0.04 

1.51 

m 

II 

0.01 

0.38 

-0.08 

-3.03 

-0.07 

-2.65 

m 

1.11 

42.01 

-0.23 

-8.70 

0.88 

33.30 

Region 

IV 

2.30 

87.04 

-0.52 

-19.68 

1.78 

67.36 

Region 

V 

-0.02 

-0.64 

-0.11 

-4.16 

-0.13 

-4.81 

Table  3.  Oil  Flow  Rates  up  the  Liner  (SAE-30  at  1500  rpm) 


mmmam 

PWR  EXCHANGE 

GAS  EXCHANGE 

TOTAL  (pwr+gas) 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

i 

-0.34 

-12.87 

-0.12 

-4.54 

-0.46 

-17.41 

H 

ii 

0.08 

3.03 

-0.21 

-7.95 

-0.13 

-4.92 

H 

in 

0.74 

28.01 

-0.23 

-8.7 

0.51 

19.3 

11111181 

rv 

1.43 

54.12 

-0.45 

-17.03 

0.98 

37.09 

V 

0.32 

12.11 

-0.19 

-7.19 

0.13 

4.92 

Table  4  Oil  Flow  Rates  up  the  Liner  (SAE-30  at  3000  ipm) 
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15W-40  at  1 500rpm 

PWR  EXCHANGE 

GAS  EXCHANGE 

TOTAL  (pwr+gas) 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

i 

-0.04 

-1.51 

0.08 

3.03 

0.04 

1.51 

fiBtBBI 

ii 

0.03 

1.14 

-0.04 

-1.51 

-0.01 

-0.38 

m 

in 

0.63 

23.84 

0.05 

1.89 

0.68 

25.73 

HHHI 

IV 

1.17 

44.28 

-0.08 

-3.03 

1.09 

41.25 

m 

V 

0.16 

6.06 

-0.01 

-0.38 

0.15 

5.68 

Table  5.  Oil  Flow  Rates  up  the  Liner  (15W-40  at  1500  rpm) 


1 5W-40  at  3000rpm 

PWR  EXCHANGE 

GAS  EXCHANGE 

TOTAL  (pwr+gas) 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow 

Vol  Diff 

Oil  Flow- 

HH 

I 

-0.22 

-8.33 

-0.16 

-6.06 

-0.38 

-14.38 

II 

0.10 

3.78 

-0.20 

-7.57 

-0.i0 

-3.78 

Region 

m 

0.46 

17.41 

-0.38 

-14.38 

0.08 

3.03 

Region 

IV 

0.91 

34.44 

-0.66 

-24.98 

0.25 

9.46 

V 

0.09 

3.41 

-0.35 

-13.25 

-0.26 

-9.84 

Table  6.  Oil  Flow  Rates  up  the  Liner  ( 1 5W-40  at  3000  rpm) 


Upon  examination  of  these  figures,  it  is  seen  that  a  majority  of  the  regions  show 
absolutely  no  correlation  to  oil  consumption.  However,  a  slight  correlation  does  exist  for 
Region  II,  tiie  top  ring,  but  only  during  the  power  exchange  strokes.  Therefore,  a  closer 
look  is  taken  at  non-dimensionalized  values  for  this  case  and  for  the  oil  consumption.  This 
is  accomplished  Py  dividing  each  rate  by  its  corresponding  1 500 
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rpm  value,  (i.e.  all  of  the  1500  rpm  values  will  equal  one)  The  results  are  illustrated  in 
figure  23  for  SAE-30  and  figure  24  for  15V/ -40,  however,  it  is  seen  that  the  oil  flow  rates 
do  not  duplicate  the  trends  of  the  oil  consumption  in  either  case. 


6.3  PROBABLE  ERRORS  IN  ANALYSIS  TECHNIQUE 

When  comparing  oil  flow  rates  to  oil  consumptions,  if  in  fact  film  thickness  in  this  case 
can  even  be  related  to  oil  consumption,  only  trends  can  be  looked  for.  An  exact  matching 
of  values  would  be  virtually  impossible  for  the  following  reasons: 

a)  Oil  consumption  measurements  vary  within  an  engine,  therefore  they  will 
vary'  even  more  between  separate,  but  identical,  engines. 

b)  The  oil  film  data  is  obtained  whiie  using  new,  and  thus  more  volatile, 
lubricating  oil.  (Refer  back  to  figure  9) 

c)  As  previously  mentioned,  the  calibration  coefficients  for  each  film  trace  are 
determined  by  way  of  etch  marks  on  the  skirt  [1].  However,  there  may 
exist  a  temperature  gradient  between  the  skirt  and  the  rings  and  lands. 

This  gradient  would  result  in  different  calibrations  for  these  latter 
regions  since  the  fluorescence  of  the  oil  is  inversely  proportional  to  its 
temperature  [6J. 

d)  The  oil  flow  relief  holes,  described  in  sect. 3. 4,  is  probably  what  attributes 
to  the  absence  of  correlation  between  oil  consumption  and  film  thickness. 
Figure  25  illustrates  that,  over  the  length  of  the  piston,  the  power  strokes 
tend  to  have  more  oil  on  the  way  up  then  on  the  way  down  and  vice  versa 
for  the  gas  exchange  strokes.  These  results  imply  that: 

1 )  after  the  oil  control  ring  passes  the  window  on  the  compression 
stroke  and  before  it  passes  the  window  on  the  expansion  stroke,  oil 
is  returned  to  the  sump  via  the  flow  holes. 

2)  after  the  oil  control  ring  passes  the  window  on  the  exhaust 
stroke  and  before  it  passes  the  window  on  the  intake  stroke,  oil 
is  supplied  to  the  sump  via  the  flow  holes. 
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6.4  MISCELLANEOUS  ANALYSIS 

The  temperature  of  the  upper  portion  of  the  liner  during  operation  is  plotted  in  figure  26 
for  the  RPMs  shown.  These  temperatures  are  then  used  to  determine  the  high  and  low 
shear  viscosities  (figure  27)  and  the  surface  tension  (figure  28)  of  the  oils  in 
use  [7],  Lastly,  the  oil  consumptions  versus  viscosities  are  plotted  in  figure  29. 
Unfortunately,  no  trends  or  relationships  between  the  two  are  observed. 


The  primary  conclusions  and  recommendations  are  as  follows: 

1 )  The  crown  land  is  virtually  dry  and  therefore  does  not  contribute  to  oil 
consumption. 

2)  The  top  ring  is  the  only  region  whose  oil  film  thickness  characteristics  remotely 
resemble  oil  consumption  However,  before  any  conclusions  can  be  formulated,  more 
comparisons  need  to  be  made  between  this  region's  flow  rate  and  oil  consumption. 

? )  The  concept  of  oil  being  "blown  ’  around  the  top  ring  and  into  the  combustion 
chamber  is  still  a  viable  solution,  however,  he  top  ring  region  might  be  an  additional 
driving  force  behind  oil  consumption. 

4)  The  oil  flow  relief  holes  provide  a  significant  unknown  which  might  make  this 
type  of  flow  rate  analysis  for  this  particular  set-up  invalid.  A  recommended  set-up  would 
consist  of  obtaining  film  thickness  measurements  through  a  window  placed  such  that  the  oil 
control  ring  never  passes  over  it 
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APPENDIX  A  DERIVATION  OF  OIL  CONSUMPTION 
CALCULATION  131 

The  following  derivation  is  based  on  the  composition  of  the  combustion  water.  This 
water  comes  from  the  water  in  the  inlet  air  and  from  its  formation  via  the  combustion 
process  of  the  fuel  and  of  the  oil.  If  radioactive  hydrogen  is  only  in  the  oil,  then: 

disintegrations  per  minute  (dpm)  =  Lube  oil  consumption  dpm  (a.  1 ) 

Wt  of  sample  of  combustion  H2O  H2O  from  fuel  comb.  +  H2O 

from  oil  comb.  +  H2O  from  air 

The  following  symbols  are  used: 

R0  =  Wt  rate  of  oil  consumption,  g/hr 
Rf  =  Wt  rate  of  fuel  consumption,  g/hr 
Ra  =  Wt  rate  of  air  consumption,  g/hr 

K()  =  Water  formation  constant  of  oil,  g  water  formed/g  oil  burned 
Kf  =  Water  fonnation  constant  of  fuel,  g  water  formed/g  fuel  burned 
hs  =  Specific  humidity,  g  H20/g  air 

SAov  =  Specific  Activity  for  H2O  sample,  dpm/g  of  sample 
SA0  =  Specific  Activity  for  lube  oil,  dpm/g 
t  =  time,  hr 

These  values,  when  substituted  into  equation,  a.l.  result  in: 

SACW  =  R0t<SA0)/[R0f  K()  +  Rfl  Kf  +  Ralh^l 
Subsequent  grouping  of  the  R0  tenns  results  in: 

R0  =  (Rfkf  +  R;,hs][(S\0'SAc„)-K0] 
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The  water  formation  constants,  k0  and  Kf.  are  determined  based  on  their  hydrogen  to 
carbon  ratios  (H/C).  (See  Appendix  E) 

The  H/C  ratios  for  the  Pennzoil  oils  are  obtained  via  Dr.  Tayeb  Benchaita  of  Pennzoil 
Products.  The  H/C  for  the  fuel  is  based  on  a  standard  light  diesel  fuel.  Actual  values  and 
sample  calculations  are  included  in  Appendix  E. 
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APPENDIX  B  DETAILED  EQUIPMENT  DESCRIPTION 

B.l  KUBOTA  EA300N  CHARACTERISTICS 


The  Kubota  EA300N's  detailed  geometry  is  presented  in  the  following  table. 


Make 

Kubota 

Model 

EA300N 

Type 

Horizontal,  4-stroke,  IDI  Diesel 

Number  of  cylinders 

1  (one) 

Bore  and  Stroke 

75  x  70  mm  (2.95  x  2.76  in.) 

Displacement 

0.309  liters  (18.86  cu.  in.  ) 

Compression  Ratio 

23:1 

Cooling  System 

Water  cooled  -  natural  convection 

Lubrication  System 

Trochoidal  pump  (no  oil  filter) 

Rated  Brake  Horsepower 

4.48  kW  (6hp)  <S>  3000  RPM 

Maximum  Torque 

15.2  Nm  (11.2  ft-lbf) 

Typical  Application 

Remote  Power  Generation 

Table  B.l  Engine  Characteristics 


A  unique  characteristic  of  this  engine  is  that  it  does  not  have  a  paper  filter.  Instead,  it 
has  a  strainer  screen  and  a  magnet  to  keep  large  pieces  of  debris  and  metal  shavings  from 
entering  the  oil  pump.  Unfortunately,  this  allows  the  dirt  and  grit  that  is  introduced  into  the 
oil  during  operation  to  remain  in  the  oil.  Therefore,  the  engine  oil  becomes  black  very 
quickly.  This  "blackness"  requires  the  oil  samples  to  be  diluted  and  centrifuged  (see 
Appendix  D)  in  order  for  the  scintillation  counting  technique  to  be  of  use. 
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I5.Z  KAUKJAlilVfc,  (ML  161 

The  tritiated  oil  received  has  been  through  a  reduction  process  using  tritium  gas.  This 
gas.  plus  one  ml  of  the  neutral  base  oil,  was  mixed  with  forty  mg  of  a  platinum  black 
catalyst  and  stirred  overnight.  Afterwards,  any  labiles  were  removed  using  50%  benzene 
methanol.  The  final  result  is  the  addition  of  tritium  atoms  to  the  hydrocarbon  molecules  in 
the  oil.  This  process  is  referred  to  as  either  a  hydrogenation  or  tritiation  process. 

B.3  ENGINE  INSTRUMENTATION 

The  following  is  a  listing  of  the  instrumentation  used  while  collecting  data  for  these 
experiments.  (Refer  back  to  Fig.  1 ) 


COMPONENT: 

Make,  Model  #: 

Liner  Thermocouples  (3) 

Omega,  J-type 

Liner  Thermocouple  Readout 

Omega,  Model  #2 1 68 A 

Other  Thermocouples 

Omega,  K-type 

Thermocouple  Readout 

Omega,  Model  #199 

Hall -effect  Speed  Sensor 

Minarik,  Visi-Tach 

Load  Cell 

Eaton,  Model  #3169 

Laminar  Air  Flow  Element 

Meriam  Instr.,  Model  50MW20-2 

Liquid  Scintillation  Counter 

Packard,  2000  CA  TRI-CARB 

Table  B.2  Measuring  Equipment  and  Instrumentation 
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APPENDIX  C  RADIOTRACER  PROCEDURES 

C.l  MIXING  AND  "BREAKINC.-IN"  THE  OIL 

Non-radioactive  oil  should  be  mixed  with  radioactive  oil  until  the  concentration  level  is 
approximately  six  microcuries  per  gram  (6  ^lCu/g).  If  a  completely  new  batch  of  oil  is 
used,  the  engine  should  be  operated  for  around  5  hours  before  any  data  samples  are 
collected.  This  break-in  period  allows  any  light  ends  to  be  burned  off.  (Refer  back  to 
figure  9) 

C.2  SAMPLING  PROCESS 

Before  samples  are  taken,  operate  the  engine  continuously  at  the  given  load  and  speed  to 
allow  temperatures  to  stabilize.  Between  speed  changes,  allow  approximately  15  minute 
for  stabilization.  After  stabilization,  run  the  vacuum  pump  for  a  short  period  of  time  in 
order  to  draw  exhaust  gas  through  the  sampling  system.  This  is  done  ( 1 )  to  ensure  that  the 
engine  and  sampling  system  are  making  good  water  samples  and  (2)  to  ensure  that  any 
water  residue  from  the  last  data  run  is  emptied  from  the  condenser  coils.  The  following  is  a 
step-by-step  description  of  the  actual  sampling  process. 

1 )  Replace  water  collection  bottle  with  a  clean  one 

2)  Replace  exhaust  filter  with  a  new  one 

3)  Turn  on  the  vacuum  pump  to  initiate  the  sampling  process 

4)  Measure  and  record  the  volume  flow  rate  of  fuel  using  the  installed  burette 
and  a  timer 

5 )  Record  the  pressure  difference  across  the  laminar  flow  element  using  the 
installed  inclined  manometer 

6)  Record  operating  temperatures 

7)  After  approximately  ten  minutes,  secure  from  sampling  by  turning  off  the 


vacuum  pump.  This  time  limit  allows  for: 
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a)  collection  of  a  sufficient  amount  of  condensed  water  and 

b)  a  sufficient  number  of  piston  ring  revolutions  and  thus  a  better 
average  oil  consumption 

8)  Remove  both  the  filter  and  the  water  collecting  bottle  and  replace  with  clean 
ones  (see  Appendix  D  for  preparation) 

9)  Stabilize  at  the  next  load  and  speed  and  repeat  from  step  3  on. 

10)  After  running  at  the  desired  load  and  speeds,  secure  the  engine  and 

drain  approximately  50  ml  of  lubricating  oil  followed  by  draining  one  ml  of  oil 
into  a  separate  test  tube  (see  Appendix  D  for  preparation) 

1 1 )  Return  the  50  ml  of  drained  oil  to  the  sump 

12)  Record  the  barometric  pressure  and  relative  humidity 

All  data  is  recorded  on  available  data  sheets  in  the  units  indicated  on  the  sheet.  Using 
these  units  allows  the  values  to  be  inputted  directly  into  the  present  spreadsheets  (Appendix 
E).  Copies  of  the  data  sheets  in  use  are  included  as  the  next  two  pages. 
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_ EXPERIMENTAL  ...DATA  SHEET 

UNIFORM  INPUT: 

DATE  - 


Oil  brand 


Oil  Activity  - 

(dpm/.5  ml  kero/oil  soln ) 


Barometric  pressure  - 

(mm  mercury) 


:er.; 


■  a  t  u  r  e 


coe.  t 


Relative  Humidity 
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Load - 

( lbs  ) 

Speed  - 

( rpm  ) 

F  '.1  e  1  rate - 

( m 1 / sec ) 

Air  press,  diff. 
across  element  - 
(inches  wtr) 

amoun t  of  water 

collected  - 

(ml  ) 

activ  of  soot 

sample  - 

i  dpm ) 

activ  cone,  of 

wtr  sample  - 

(dpm /ml )  or  ( dpm/g ' 

TEMPERATURES: 

exhaust  - 

air  inlet  - 

oil  - 

coolant  - 


upper  - 

middle  - 

lower  - 


Load - 

(lbs) 

Speed  - 

(  rpm  ) 

Fuel  rate  - 

(ml /sec ) 

Air  press .  dif f  . 
across  element  - 
(inches  wtr ) 

amount  of  water 
collected  - 

/  -n  J  \ 

t  III  1  ,■ 

activ  of  soot 

sample  - 

(dpm  ) 

activ  cone .  of 

wtr  sample  - 

(dpm /ml)  or  (dpm/g) 

TEMPERATURES: 

exhaust 

air  inlet  - 

oil  - 

coolant  - 

liner: 

upper  - 

middle  - 

lower  - 
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APPENDIX  D  SAMPLE  PREPARATION  [71 

D.l  WATER  SAMPLES 

1 )  Filter  the  condensed  water  into  a  clean  test  tube 

2)  Add  and  mix  a  sufficient  amount  of  ammonium  hydroxide  solution  to  the 
sample  to  bring  the  PH  up  to  7  (neutral).  [7,8] 

3)  Centrifuge  the  sample  for  approximately  ten  minutes 

4)  Filter  again  if  necessary 

5)  Pipette  10  ml  of  fluor  solution  in:o  a  counting  bottle 

6)  Pipette  1  ml  of  the  water  sample  into  the  counting  bottle  (note  -  more  or 
less  water  may  be  counted  but  the  spreadsheet  is  set  up  for  1  ml  —see 
Appendix  E) 

7)  Replace  top,  mix,  and  number  the  sample 

P.2  SOOT  SAMPLES 

1 )  Place  the  filter  in  a  tube  containing  approximately  30  ml  of  fluor  solution  and 
shake  vigorously  to  mix  well 

2)  Allow  two  hours  for  the  fluor  solution  to  dissolve  the  soluble  hydrocarbons 
and  then  shake  again 

3)  Filter  the  solution  into  a  clean  test  tube.  Repeat  if  necessary 

4)  Pipette  10  ml  of  the  fluor  solution  into  three  counting  bottles  each 

5)  Divide  the  filtered  solution  among  the  bottles 

6)  Replace  the  tops,  mix,  and  number  the  samples 

D.3  OIL  SAMPLES 

1 )  Pipette  10  ml  of  kerosene  into  a  clean  test  tube 

2)  Pipette  .2  ml  of  the  oil  sample  into  the  kerosene.  Subsequently,  flush 
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the  pipette  with  the  oil/kerosene  solution  until  all  visible  oil 

residue  is 

removed  from  the  pipette.  Place  pipette  into  the  solution 

3)  Mix  well  and  centrifuge  for  approximately  ten  minutes 

4)  Pipette  10  ml  of  the  fluor  solution  into  a  counting  bottle 

5)  Pipette  .5  ml  of  the  oil/kerosene  solution  into  the  counting  bottle  (note  - 
as  before,  more  or  less  may  be  counted  but  the  spreadsheet  is  set  up  for  .5 
ml  —  see  Appendix  E) 

6)  Replace  top,  mix.  and  number  the  sample 

All  samples  are  taker,  to  the  Radiation  Protection  Office  so  that  their  activity  levels  may 
be  counted  using  a  liquid  scintillation  counter. 
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APPENDIX  E  FORMULAS  AND  EXAMPLE  CALCULATIONS 

E.l  OIL  CONSUMPTION 

E.Ia  MEASURED  CONSTANTS 

pf  =  fuel  density,  g/ml  =  0.838 

po  =  oil  density,  g/ml  =  0.826  (15V/ -40)  or  0.841  (SAE-30) 

(H/C)f  =  hydrogen  to  carbon  ratio  of  fuel  =  1.8  (light  diesel) 

(H/C)0  =  hydrogen  to  carbon  ratio  of  oil  =  1.88  (15W-40)  or  1.78  (SAE-30) 
LFE  =  Laminar  flow  element  rating.  CFM/in  water 

E.lb  CALCULATED  CONSTANTS 

Kf  =  water  fonnation  constant  of  fuel 
=  [g  H20/gH]x[g  H/gfuel] 

=  [(2  +  16)/2]  x  f( H/C )f/(  1 2  +  (H/C)f)] 

=  1.174 

K.0  =  water  formation  constant  of  oil 
=  1.219  (15W-40)  or  1.163  (SAE-30) 


E.Ic  MEASURED  VARIABLES 
Vf  =  volume  flow  rate  of  the  fuel,  ml/sec 

dP  =  differential  pressure  across  laminar  air  flow  element,  inches  water 

\Y  =  amount  of  water  condensed,  ml 

SAov'  =  specific  activity  of  condensed  water,  dpm/ml 

As  =  total  activity  of  soot  sample 

SAo'  =  specific  activity  of  oil/kerosene  mixture,  dpm/(.5  ml  of  oil/kero  solution) 
atm  =  barometric  pressure,  mm  of  mercury 
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Ta  =  air  temperature,  °C 
H  =  relative  humidity,  % 

E.ld  CALCULATED  VARIABLES 

hs  =  specific  humidity,  g  H20/g  air 

=  determined  from  psychrometric  chart  using  H  and  Ta 
Rf  =  wt  rate  of  fuel,  g/sec 
=  Vfxpf 

(Xc  =  viscosity  coefficient 

=  1  -  [.0025714  x  (Ta  -  20)]  (note  —  accurate  for  15  <  Ta  <  35) 
pa  =  density  of  air,  g/m? 

=  (atm)  x  (133.32  Pa/mm  Hg)  x  (28.962  g/mol)  x  (mol  °K/8.3143  Nm)  x 
[1/(273  +  Ta)°K]  (note  —  ideal  gas  law) 

Ra  =  wt  rate  of  air,  g/sec 

=  (LFE)  x  (dP)  x  (.02832  nvVft3)  x  (min/60  sec)  x  [(atm)/29.92  in  Hg]  x 

(in/25.4  mm)  x  [530  °R/[460  +  (1.8Ta  +  32)]°R)  x^-xpj,  [8] 
SAC\y  =  total  specific  activity  of  exhaust  sample,  upm/g 
=  SAcw*  +  (Ajj/W  ) 

SAo  =  total  oil  consumption,  g/hr 

=  [RfKf  +  Rahs]/[(,SACW/S A0)  -  K0]  x  (3600  sec/hr) 

E.2  OIL  FLOW  RATES  (FROM  FILM  THICKNESS  TRACES) 

RPM  =  revolutions  per  minute 

V()p  =  average  oil  volume  difference,  for  a  given  region,  between  the 
power  exchange  strokes,  mm-  /cycle 
=  calculated  using  an  existing  computer  program  which  was  modified 


for  the  Kubota  (Appendix  F) 
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V'oji  =  average  oil  volume  difference,  for  a  given  region,  between  the 
gas  exchange  strokes,  mm-Vcycle 

=  calculated  using  an  existing  computer  program  which  was  modified 
for  the  Kubota  (Appendix  F) 

V0  =  total  avg.  volume  difference  between  the  up  and  down  strokes,  mm- 

=  V op  +  V ()g 

Hop  =  rate  of  oil  flow  up  the  liner  based  on  film  thicknesses  during  the 
power  exchange  strokes,  g/hr 
=  (V0p)  x  (cnfVlOOO  mm^)  x  (RI»M)  x  (60  min/hr)  x 
(cycle/2  revolutions)  x  (p0) 

Rog  =  rate  of  oil  flow  up  the  liner  based  on  film  thicknesses  during  the 
gas  exchange  strokes,  g/hr 

=  (Yog)  x  (cm3/ 1000  nun3)  x  (Rpjvj)  x  (60  min/lir)  x 
(cycle/2  revolutions)  x  (p0) 

Ro'=  overall  rate  of  oil  flow  up  the  liner  based  on  film  thicknesses,  g/hr 
=  Rop  +R«g 

All  of  the  preceding  formulas  are  incorporated  into  LOTUS  123  v2.0  spreadsheets. 
Examples  of  these  spreadsheets  are  included  as  the  following  six  pages.  The  first  two 
pages  are  for  oil  consumption  calculations  and  the  last  four  are  for  oil  rates  based  on  the 


film  thickness  data. 
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E.3  LOTUS  123  v2.01  SPREADSHEETS 


OIL  CONSUMPTION  CALCULATIONS: 

UNIFORM 


INPUT  CALCULATIONS 


oil  H/C  ratio  - 

1.78 

/ 

Oil  formation  constant  - 

--1. 162554 

( SAE-30  =  1.78) 

/ 

( 15W-40  =  1.88) 

/ 

/ 

fuel  H/C  ratio  - 

1.8 

/ 

Fuel  formation  constant 

--1 . 173913 

/ 

activ  of  oil  smpl  - 

57052 

/ 

activ  cone  oil  - 

— 6919505 . 

(dpm/.5  ml  soln) 

/ 

( dpm/g ) 

/ 

oil  density  - 

0.841 

/ 

(SAE-30=  .841) 

/ 

( 15W-40= . 826) 

/ 

/ 

Barometric  Pressure  - 

777 

/ 

(mm  mercury) 

/ 

/ 

humidity  coeff  - 

0.01 

/ 

fib  wtr/lb  air) 

/ 

/ 

air  temperature  - 

28 

/ 

( Celsius ) 

/ 

1 


UNIQUE 


INPUT 

Load  (lbs)  - 

Speed  (rpm) - 

fuel  rate  - 

(ml/sec ) 

air  pressure  diff  - 
(inches  water) 


activ  of  soot  smpl - 

(dpm/1  ml  soln) 

activ  cone  wtr  smpi---- 
(dpm/ml)  or  (cipm.'cr) 

amount  wtr  collected  -~ 
( ml ) 


7.6 

1500 

**»  n,  ■*>  O 
•  ^  V  w 

1.27  5 


6  3  ? 


n: 


/  CALCULATIONS 

/ 

/  TOTAL  oil  consumption  - 1.425019 

/  (g/hr) 

/ 

/ 


/  fuel  rate - 0.170734 

/  (g/sec) 

/ 

/  air  rate - 3.529377 

/  (g/sec) 

/ 

/  equivalence  ratio - 0. 693130 

/  (F/A)S  =  0.0693 
/ 

/  total  activ  of  soot  smpl  -  632 

/  ( dpm ) 

/ 


/  total  cone,  of  sample  - 11593.93 

/  (dpm/ml'i  or  (dpn/g) 

/ 

t 

I 

/ 

/ 


/  oil  cor.sv.mp  (burned)  - 1.39313 

/  oil  con sump  (unburn)  - 0.02673 

/  Total - 1.424916 


.u 
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OIL  FLOWS  UP  THE  LINER: 


PISTON 

REGIONS 

AVG  VOL 

OIL  FLOW 

AVG  VOL 

OIL  FLOW 

SAE-30 

PWR  STRK 
( cu  mm) 

PWR  STRK 
(g/hr) 

GS  XCHNG 
{cu  mm) 

GS  XCHNG 
(g/hr) 

1500 . 00 

I 

II 

1/2  crwn  Ind 
top  ring 

(X2)  0.04 

0 . 01 

1.51 

0.38 

0.00 

-0.08 

0.00 

-3.03 

III 

1  ring  width 

below  1.11 

42.01 

-0 .23 

-8.70 

IV 

2nd  land 

2.30 

87.04 

-0.52 

-19 . 68 

V 

2nd  ring 

-0.02 

-0 .64 

-0.11 

-4 . 16 

PISTON 

AVG  VOL 

OIL  FLOW 

AVG  VOL 

OIL  FLOW 

REGIONS 

PWR  STRK 

PWR  STRK 

GS  XCHNG 

GS  XCHNG 

SAE-30 

{cu  mm) 

( g/hr) 

(cu  mm) 

(g/hr) 

3000 . 00 

I 

1/2  crwn  lnd 

{ x2 )  -0.34 

-12.87 

-0.12 

-4.54 

II 

top  ring 

0 . 08 

3 . 03 

-0.21 

-7.95 

III 

1  ring  width 

below  0.74 

28.01 

-0.23 

-8.70 

IV 

2nd  land 

1  .43 

54  .12 

-0.45 

-17 . 03 

V 

2nd  ring 

0.32 

12 .11 

-0 .19 

-7  .  19 

PISTON 

AVG  VOL 

OIL  FLOW 

AVG  VOL 

OIL  FLOW 

REGIONS 

PWR  STRK 

PWR  STRK 

GS  XCHNG 

GS  XCHNG 

15W-40 

{ cu  mm ) 

( g/hr) 

(cu  mm) 

( g/hr) 

1500 . 00 

T 

1/2  crwn  lnd 

( x2 )  -0.04 

-1.51 

0.08 

3.03 

II 

top  ring 

0.03 

1-14 

-0.04 

-1.51 

III 

1  ring  width 

below  0.63 

23.84 

0.05 

1.89 

IV 

2nd  land 

1  .17 

44.28 

-0.08 

-3 . 03 

V 

2nd  ring 

0  .  16 

6.06 

-0 .01 

-0 .38 

PISTON 

AVG  VOL 

OIL  FLOW 

AVG  VOL 

OIL  FLOW 

REGIONS 

PWR  STRK 

PWR  STRK 

GS  XCHNG 

GS  XCHNG 

15W-40 

( cu  mm ) 

( g/hr) 

(cu  mm) 

(g/hr) 

3000 . 00 

I 

1/2  crwn  lnd 

( X2 )  -0.22 

-8.33 

-0 .16 

-6 . 06 

II 

top  ring 

0  .  10 

3.78 

-0.20 

-7  .  57 

III 

1  ring  width 

below  0.46 

17.41 

-0.38 

-14.33 

IV 

2nd  land 

0 . 91 

34 . 44 

-0 .66 

-24  .98 

V 

2nd  ring 

0 . 09 

3.41 

-0.35 

-13 . 25 
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OIL  FLOWS 


SAE-30 
1500 . 00 


SAE-30 
3000 . 00 


15W-40 
1500 . 00 


15W-40 
3000 . 00 


UP  THE  LINER: 


PISTON 

AVG  VOL 

OIL  FLOW 

REGIONS 

ALL  STRK 
(cu  mm) 

ALL  STRK 
(g/hr) 

I 

1/2  crwn  lnd 

(x2)  0.04 

1.51 

II 

top  ring 

-0.07 

-2.65 

III 

1  ring  width 

below  0.88 

33.30 

IV 

2nd  land 

1.78 

67.36 

V 

2nd  ring 

-0 . 13 

-4.81 

PISTON 

AVG  VOL 

OIL  FLOW 

REGIONS 

ALL  STRK 
( cu  mm) 

ALL  STRK 
(g/hr) 

I 

1/2  crwn  lnd 

( x2 )  -0.46 

-17.41 

II 

top  ring 

-0 . 13 

-4 . 92 

III 

1  ring  width 

below  0.51 

19.30 

IV 

2nd  land 

0.98 

37 . 09 

V 

2nd  ring 

0 . 13 

4.92 

PISTON 

AVG  VOL 

OIL  FLOW 

REGIONS 

ALL  STRK 
(CU  mm) 

ALL  STRK 
(g/hr) 

I 

1/2  crwn  lnd 

(x2 )  0.04 

1.51 

II 

top  ring 

-0.01 

-0.38 

III 

1  ring  width 

below  0.68 

25.73 

IV 

2nd  land 

1 . 09 

41 . 25 

V 

2nd  ring 

0.15 

5.68 

PISTON 

AVG  VOL 

OIL  FLOW 

REGIONS 

ALL  STRK 
( cu  mm) 

ALL  STRK 
(g/hr) 

I 

1/2  crwn  lnd 

( x2 )  -0.38 

-14.38 

II 

top  ring 

-0 . 10 

-3.78 

III 

1  ring  width 

below  0.08 

3.03 

IV 

2nd  land 

0.25 

9.46 

V 

2nd  ring 

-0.26 

-9 . 84 
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APPENDIX  F  DETERMINATION  OF  VOLUME  DIFFERENCES 

F.l  DETERMINATION  OF  REGIONS 

The  program  utilized  uses  a  Crank -Nicolson  technique  to  integrate  the  difference  in 
areas  under  a  film  trace  for  the  upstroke  and  downstroke  pairs  of  specific  regions. 
Subsequently,  these  area  differences  are  swept  around  the  circumference  of  the  bore  to 
determine  the  volume  differences. 

Each  inputted  region  is  calculated  based  on  using  the  "zero-point"  crank  angle  as  a 
reference  position.  This  position,  initially  detennined  to  equal  -81.6°,  is  defined  as  the 
angle  at  which  the  top  of  the  piston  passes  the  quartz  window  on  an  upstroke.  The 
following  relation  is  then  used  to  determine  the  distance,  s,  in  cm,  from  the  crank  axis  to 
the  piston  pin  axis  at  this  point: 

s  =  a  x  cos(0)  +  (l2  -  (a2  x  sin2!©))1/2 

where:  1  =  connecting  rod  length  =  1 1 .01  cm 
a  =  crank  radius  =  3.5  cm 
0  =  crank  angle 

Next,  a  new  "zero-point"  crank  angle  is  obtained  by  fitting  a  model  of  the  piston  to  a 
printout  of  the  film  trace.  Since  piston  geometry  is  known,  this  step  allows  for  an  estimate 
in  the  error  of  the  original  "zero-point".  This  error  distance  is  subtracted  (or  added)  to  the 
original  s  and  a  new  "zero-poinf’crank  angle  is  detennined  by  iterating.  Specific  regions 
can  then  be  determined  in  a  similar  manner. 

For  example: 

a)  the  initial  "zero-point"  is  -81.6°.  Therefore,  s0  =  10.96  cm 

b)  it  is  determined  from  the  trace  that  the  s  is  0.05  cm  off.  Therefore, 
the  new  s0  =  1 0.96  -  0.05  =  10.91  cm 

c)  by  iterating,  the  new  "zero-point"  is  determined  to  equal  -82.5° 

d)  the  crown  land  is  0.95  cm  long 

e)  therefore,  the  new  s  =  sj  =  sp  +  d  =  10.94  +  0.95  =  1 1.89  cm 
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d)  now  iterate: 

0  =  -67.7°  implies  Si  =  1 1 .85 
0  =  -67.9°  implies  =  1 1 .84 
0  =  -67.5°  implies  Sj  =  1 1 .86  (correct  one) 

e)  the  region  relating  to  the  crown  land  would  then  be  -82.5°  to  -67.5° 

f)  the  top  ring  region  is  determined  in  a  similar  manner  except  now 
s0  =  11.91  cm 

Regions  n,  IV,  and  V  consist  of  the  top  ring,  second  land,  and  second  ring, 
respectively.  They  are  used  to  reflect  the  film  behavior  of  the  upper  portions  of  the  piston. 
The  crown  land  is  not  included  based  on  the  reasoning  presented  in  section  6.2.  Region  I 
is  determined  based  on  the  film  on  the  liner  as  three  centimeters  of  the  crown  land  pass  the 
window.  The  assumption  is  that  the  thickness  on  the  liner  at  this  point  is  the  same  up  to 
TRR  (top  ring  reversal).  Finally,  region  IE,  one  ring  width  below  the  top  ring,  is  used  in 
an  attempt  to  substantiate  the  claim  that  oil  is  drawn  up  and  around  the  top  ring  during 
operation. 


ooooooooooooooooooooooonoonoononooooooooooooo 


46 


F.2  "MASSBAL . FOR”  LISTING 

C 

PROGRAM  MASSBAL 
C 

C  •««  Program  to  compute  up  /  down  stroke  oil  volume  di'ferences  in 
C  distinct  specified  regions  along  the  piston  for  multiple  Sloan 

C  data  files  and  write  them  to  an  output  file  that  can  be 

C  converted  to  a  bar  graph  format. 

C 

C  This  program  requires  as  input  any  number  of  standard  Sloan  data 
C  files  containing  an  array  of  position  vs  film  thickness.  The  program 

C  will  average  the  volume  of  oil  in  any  number  of  user  specified 

C  regions  along  the  piston, 

c 

c  The  program  has  been  adapted  to  analyze  either  McElwee 

c  [where  IRQA7E(3)-39 ]  or  Bliven  flor.»  data  sets.  Commands 

c  specific  to  Lux  [where  IRDATE(3) -£8]  data  sets  have  been 

c  neutralized  by  comment  lines  which  are  specified  by  "cc". 

c  [M.  Bliven,  14  March  1990] 


•••  Progrcm  structure: 

Main  Level  1st  Sublevel  2nd  Sublevel  3rd  Sublevel 

MASSBAL  FSTOCTPT 

INTECRATEDF 

AVERAGE 

...  I/O 


Uni  t 

3: 

INFILE  Sloan 

data  file 

Un  i  t 

5: 

Terminal 

I 

Un  i  t 

6: 

Terminal 

0 

UNIT 

10: 

OUTFILE. 

BDC 

LOCATIONS.  RPM 

• • •  impor t an t  variables 

AVALS  Array  of  average  volume  for  like  revolutions 

CALOC  Crank  angle  location 

NPMAX  parameter,  max  number  of  points  to  plot 

CALFC7R  calibration  factor  for  individual  dato  set 

CACFFS7  correction  factor  for  shaft  encoder  offset 

FILM(2.NFVAX)  Array  of  position  vs  film  thickness 
GEMNVAL  Mean  value  of  oil  volume  difference  for  gas 

exchange  strokes 

GESO  Standard  deviation  for  oil  volume  difference 

for  gas  exchange  strokes 
INFILEV  Input  file  vector 

INTRFM  Value  for  engine  speed  calculated  from  time 

difference  between  BDC  pulses 
LWR8(10)  Lower  bound  for  region  i 

MNVAL  Mean  value  of  oil  volume  for  a  region  for  like 

revolutions  from  the  progrom  AVERAGE 
NBDCLOC  Botton  dead  center  location  in  crank  angle 

NS7RCKE  Number  of  strokes  recorded  in  data  file 

NUMFIL  Number  of  Sloan  files  to  process 

NUMREG  Number  of  regions  to  divide  data  into 

OFFST  Instrument  offset  factor  to  initialize  zero 

f i Im  thickness 

FYiVNVAL  Mean  value  of  oil  volume  difference  between 

powered  strokes 

PWSD  Standard  deviation  or  oil  volume  difference 

between  powered  strokes 


o  o  o  o  o  o  o 
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REGVAL  Array  containing  piston  regions  with  the 

average  value  of  oil  volume  difference  for 
like  revo I u t i ons 

SIGVAl  Standard  deviation  of  oil  volume  for  like 

revolution  from  the  program  AVERAGE 


PARAMETER  NPVAX  -  190000  I  max  f  of  points/ch 

PARAMETER  NPCMAX  -  60000  1  max  f  of  points/eycle 

PARAMETER  NCMAX  -  250  1  max#  of  cycles 

PARAMETER  NREGMAX  -  15 
C 

INCLUDE  1 SLNCOM. FOR'  I  contains  IMPLICIT  INTEGER. A 

C 

INTECER.2  IHED1 (256) . IHED2(256) 

EQUIVALENCE  ( IHED1 , NRUN) , ( IHED2 . NSLOT) 

C 

INTEGERS  IDATA(NPUAX) 

C 

INTEGER  NBDCLOC( NCMAX) .NCHOICE.  NSTROKE(NCMAX) , 

$  NUMREV , PCH , ICH ,Q , ICH1 , CAL0C(365) , NLVREG 

REAL  P I LM ( 2 , NPCMAX ) ,  AVALS (2 , NCMAX) ,  VALUE(2) , 

$  XOLD.XNEW,  INTRPM(NCMAX) ,  CALFCTR , ROLOFF , CAPPT , 

S  XSCL(-i)  .  INTAR(  10000) ,OUTAR( 10000) .DUMCT. 

S  SUM.  VAL.  CFFST.LMVAL.  VMPOS .  SIGVAL.  SIGPCS . DUMP LOT. 

$  LSTROKE.  CKRAD.  CROL.  BORE.  CAOFFST.  LWRB(NREGUAX) . 

$  UPRB(NREGMAX) .  PWLWVAL( NREGMAX) ,  PWSO(NREGMAX) , 

$  CEMNVAL( NREGMAX) .  GESD(NREGMAX) ,  REGVAL( NREGMAX .NCMAX) 

C 

CHARACTER  INFILEV(20).A0,CUTFILE.a0,A»1 .9*1 ,C«1 . 

S  CALINT.3 . STROK£«4 , sk i p»3 
C 

C  •••  initialize: 

C 

CALINT  -  'ON' 

ICH  -  3 
C 

C  •••  sot  up  menu: 

C 

10  CONTINUE 

ISTAT  -  LIBSERASE_PAGE( 1 , 1 )  I  erase  page... 

WRITE  (5,1000)  CALINT 
READ  (5. ..ERR  -  10. END  -  999)  NCHOICE 
C 

C  ...  NCHOICE  «  1:  Sot  data  files  to  read 
C 

IF  (NCHOICE. EQ. 1)  THEN 

WRITE  (6,'(’’  How  many  f I o r . •  files?  (up  to  20):  *',$)') 
REA0(5 , • )  NUMFIL 
DO  1-1 .NUMFIL 

WRITE  (S.’C'  Enter  filename  f  ”$)’)I 

READ(5, ’ (A40) ' )  INFILEV(I) 

ENCDO 
END  IF 
C 

C  NCHOICE  -  2:  Set  piston  regions  by  crank  angle 
C 

IF  (NCHOICE. EQ. 2)  THEN 

WRITE  (6.’(''  How  many  piston  regions?  (up  to  16):  '',$)’) 
READ(5. •)  NUMREG 
DO  1-1. NUMREG 

WRITE  (6,’(’’  Enter  lower,  then  upper  bound  (CA  deg)  for 
$  region  /  ",I2,M:  ”$)’)! 


ooooooo  ooo 
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READ(5 , •)  LWRS(I).  UPR8(I) 
ENCCC 
ENOIF 


•••  NCHOICE  -  3:  Linearly  interpolate  between  crank  angles 


$ 


IF  (NCHOICE. EQ. 3)  THEN 

WRIT£(6, ‘ (1X‘ 'Linearly  interpolate  between  crank  angles? 

(y/n)  ".*)•) 

READ( 5 , ’ (A1 ) ' )  C 

IF  (C.EQ. ' Y ' .CR.C.EQ. * y * )  CALINT-'ON' 

IF  (C.EQ. ‘N‘ .CR.C.EQ. 'n’)  CALINT- 'OFF ’ 

ENOIF 


•••  NCHOICE  -  4:  Check  entered  information  ••• 

Note:  The  running  time  of  the  program  could  be  lengthy 
depending  on  the  number  of  data  files  and  regions  specified.  It  is 
wise  to  cheek  the  input  data  for  accuracy  prior  to  execution. 

IF  (NCHOICE. EQ. 4)  THEN 

ISTAT  -  LI8SERASE_PACE( 1 . 1 )  lerase  page 
DO  I  -  1 .  NOWFIL 

WRITE  ( 6  .  '  (  "  Data  File  f  ".13,"  is:  ".A40)')I. 

$  INFILEV(I) 

ENDOO 

DO  I  ■  1,  NLMtEG 

WRITE  (6,'("  Crank  angles  bounding  region  /".I2. 

$  "  are:  ".F7.2."  to"  .  F7 . 2)  *  )  I  .  LWRB(I),  UPRB(I) 

ENCDO 

WRITE  (6.'("  RTN  to  continue  ".$)') 

READ(5 . ' (A1  )  ' )  3 
ENOIF 
C 

C  •••  NCHOICE  -  5  Execute  main  progrom  and  write  output. 

C 

IF  (NCHOICE. EQ. 5)  THEN 
DO  I  -  1 .  NUUFIL 

ISTAT  -  LIBSERASE_PAGE( 1 . 1 )  lerase  page 
WRITE  (5.1001)  CALINT,  INFILEV(I) 

ILUN  -  3 

OPEN(UNIT-I LUN .NAME- INFI  LEV(  I ) .  TYPE-  ’OLD’  .ACCESS- 'DIRECT ' 

$  .  FOR)*- ‘UNFORMATTED  ’ ) 

C 

C  •••  and  read  header: 

C 

READ  ( I LUN ' 1 )  (IHEDI(L).L-I .256) 

READ  ( I LUN ' 2 )  ( IHED2(L) . L-t ,  256) 

C 

I DOM  -  1 

MPPC  -  IDUM.N8P0256 
C 

C  •••  load  up  the  private  common 
C 

C  These  values  are  user  specified  prior  to  data  collection. 


CAPPT  -  RUSER(2S) 

1  crank  angles  per  point 

LSTROKE  -  RUSER(27) 

1  engine  stroke  1 

‘cm' 

8CRE  -  RUSER(28) 

I  engine  bore  | 

cm' 

CRDL  -  RU$ER(29) 

I  connecting  rod  length  | 

.cm! 

C 

C  •••  Convert  to  mm 
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c 

BORE  -  bore*  10. 

LSTRCKE  -  LSTROKE* 1 0 . 

CRDL  -  CROL* 1 0 
CKRAO  -  LSTROKE  /2. 

r* 

C  •*•  Set  shaft  encoder  error  in  variable  CAOFFST 
C 

C  The  shaft  encoder  marks  the  crank  angle  at  BDC.  Depending  on 

C  certain  errors  the  recorded  BOC  position  is  not  necessarily  at  —180. 

C 

cl  IF  ( IRDATE(3) . EQ. 88)  THEN 

ce  CAOFFST  -  -178.5 

cc  ELSE 

if  ( i r da t e (3) . eg . 89 . and . nrun . g t . 99 . and . n run . I t . 1 06)  then 
caoffst  -—181 
else 

if  ( i rda t e(3) . eq . 90)  then 
caoffst  -  —180.926 

else 

CAOFFST  -  -180 
end  i  f 
end  i  f 

cc  ENOIF 

C 

C  •••  Look  for  old  BOC  file.  The  BDC  file  contains  the  specific 
C  calibration  factor  ar.d  the  number  of  revolutions  in  the  data  sot. 

C 

OUTFILE  -  1  1 

ENCOC£(  40  ,  '  ( 1  ' BOCLOC .  •  ‘  .  13) 1  .OUTFI  LE)NRL'N 
C 

C  •••  Open  existing  BOC  file 
C 

OPEN(UNI T- 1 0 . NAME1— DUTF I LE , TYPE- ' OLD ' .ACCESS- ' SEQUENTIAL ' . 

S  FORM-  ‘  FORMATTED '  ) 

READ(  1  0  .  •  )  CALFCTR,  NUV«EV 
00  J-1 .NCUAX 

REAO ( 1  0  ,  •  . ENO“ 1 3 , ERR— 3)  NBOCLOC(J). INTRFV( J ) . NSTRCKE( J  ) 

ENCDO 

15  CONTINUE 
C 

C  Sot  region  along  piston  to  analyze 

C  ' 

00  N  -  1  .  2 
DO  J-1  .NUk«EG 
DO  L-1  ,NUV«EV 

REGVAL(J.L)  -  0.0 
ENDOO 
ENOOO 
W  -  1 

DO  K— N .  NUVREV-1 .  2 
I FPT-NBDCLOC (K) 

C 

C  •••  Set  IFPT  and  NRPTS  Set  PCH  for  lobels 
C 

NRPTS  -  NBDCL0C(K+1 )  -  NBDCLOC(K)  +1  I  f  of  points  to  read 
IF  ( NRPTS. GT.NPUAX)  NRPTS  -  NPUAX 
C 

C  Linearly  interpolate  between  crank  angles 
C 


C 


IF  (CALINT. EQ. 'OFF' )  GOTO  21 
ICH1  -  1 


ooo  oooooo 
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CALL  FSTDGTPT ( I LUN . ICH1 , I FPT , NRPTS . IDATA) 

C 

NUMCA  -  0 

XOLD  -  CSCL(ICH1)«IDATA(1)+CBIA(ICH1) 
c2  IF  (NRUN .LT.4-7.ond. i rdate(3).eq.88)  XOLD-XOLD-5.0 

cc  IF  ^NkuN . GT . 4/ . ana . i raa t e (5) . eq . 26)  XCLD-aOLDtS . 3 

skip  —  'no' 

C 

DO  L-2.  NRPTS 

XNEX  -  CSCL( ICH1  ) • IDATA( L)+CBIA( ICH1 ) 
c3  IF  (NRUN . LT . 47 . and .  i  rdote(3).eq.88)  XNEW-XNEW-5 . 0 

cc  IF  (NRUN  . GT . 47 . and .  irdate(3).eq.38)  XNEYf-XNEVH-5 . 0 

IF  ( XN  Eitf .  GT .  2 . 5 . ANO . XOLD . LT . 2 . 5 )  THEN 
if  (copp t . eq . . 5 . and . sk i p . en . ' y es ’ )  then 
skip-  ' no  1 
go  to  20 
end  i  f 

if  (cappt.eq..5.and.skip.eq.‘no  ’)  skip-* yes’ 

NUMCA— NUMCA+1 
CAL0C( NUMCA )-L  +  IFPT 

20  END  IF 
XOLD-XNEW 

ENDDO  i  END  OF  L  LOCR  FOR  CRANK  ANGLE  LOCATION 

21  CONTINUE 
C 

C  •••  Avaraga  first  100  data  points  of  pmt  signol  to 
C  astaDlish  zero  offset.  .005  sets  this  to  2  microns 

C 

cl  NOTE:  This  section  adapted  to  analyze  data  collected 

cc  by  8LIVEN  whose  PVfT  signal  due  to  instrument  offset  and 

cc  laser  reflectances  alone  (i«.  dry  liner)  was  measured  to  be 

cc  approximately  0.008  V.  Averaging  over  first  100  data  points 

cc  is  by-passed  and  offset  value  is  set  to  0.008. 
cc 

CALL  FSTDGTPT (I LUN. ICH, IFPT. NRPTS, IDATA) 
c 

if  ( i rda t e(3) . eq . 89)  then 
SUM  -0. 

DO  L  -300.  400 

VAL  -  CSCL(ICH).IDATA(L)  +  C8 I A( ICH) 

SUM  -  SUM  +  VAL 
DLMCT  -  L-299 

ENCOO  I  END  OF  L  LOCR  FOR  OFFSET 
c 

OFFST  -  SLV/DUMCT-.016  I  zero  pt  offset  for  UcElwee 
e  I  se 

OFFST  -  .0  I  zero  pt  offset  for  Bliven 

end  i  f 


•••  Convert  the  data  and  write  them  into  plot-array 

Note:  CLKRTE  is  clock-rate  in  [kHz]  and  con  be  used  to 
convert  x— axis  data  time 


DUVPLOT-  CAOFFST 
0-1 

DO  L  -  1  .NRPTS 

Convert  time  to  crank  angle  by  speed  or  position  interpolation 
IF  (CALINT.EO. ’OFF')  THEN 

FILM(I.L)  -  CAOFFST  +  360 . »FL0AT(L-1 )/FL0AT(NRPTS-1 ) 

ELSE 
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C 

C  •••  If  Crank  angle  pulse  is  coincident  with  sample  number,  then 
C  increment  x  value  by  one  crank  angle 
C 

IF  ( ( L+I FPT ) .  EQ . CALOC(Q) )  THEN 
FILM(1.,.',  -  CAOFFST  +  FLOAT (C) 

DUMPLOT-FILM(I.L) 

Q-Q+1 

ELSE 

C 

C  *•»  Otherwise  divide  fractions  of  crank  angle  evenly  between  pulses 
C 

IF  (Q.EQ.1)  THEN 

FILM(I.L)  -  CAOFFST-rFLOAT ( L-1  )/FLOAT (CALCC(Q)-I FPT ) 

ELSE 

FI LM( t . L)-CLUPLOT+1 . 0/FLOAT ( (CALCC(Q)-CALOC(O-I ))) 

DUMP LOT- F I LM( 1 , L) 

END  IF 
END  IF 
END  I  F 

FI LM(2 . L)-(CSCL( ICH) • I0aTA(L)+C3IA( ICH)-OFFST)/CALFCTR 

ccc 

ccc  wr  i  te(S.  ■  ( • -FILW(2.L)  -  •\M0.5)‘)  film(2,l) 

ccc  write(6,’(2iU,i6,2f8.0)’)n,k,l,filiti(1,l),film(2,l) 

ccc 

ENOOO  I  ENO  OF  L  LOOP  FCR  FILM  ARRAY 
C 

C  •••  Call  integration  routine  for  each  region.  Program  INTEGRATEDF 
C  will  calculate  the  volume  of  oil  in  a  specified  region  for  one 
C  revolution. 

C 

DO  J-1.  NUMREG 

CALL  INTEGRaTEDF(FILM.LWR9(J) .UPRS(J) . VALUE. NRPTS.CHCL, 

$  CXRAO ) 

RECVAL(J.U)  -  VALUE(2)-80REO.  14159 

ccc 

ccc  wr i te(6. ’ (4i5.5«i0.4) ’ )m,n.k, j . value(1  ) . volue(2) . regval ( j ,m) 
ccc 

ENCOC  I  END  J  LOOP  FOR  REGIONS 
AVALS( 1 ,M)  -  VALUE( 1 ) 

M-a+i 

ENCCO  1  END  OF  K  LOOP  FOR  ALL  REVOLUTIONS  OF  1  TYPE 
C 

C  •  ••  Call  average  routine  for  like  revolutions  of  one  region.  Program 
C  AVERAGE  will  calculate  the  average  volume  for  the  like  revolutions. 
C 

DO  J-1 .NUMREG 

DO  L  -  1  .  NUMREV 
AVALS(2,L)  -  0. 

ENOOO 
DO  L-1  ,U 

AVALS(2.L)  -  REGVAL(J.L) 

ENODO 

CALL  AVERAGE ( AVALS, fcWVAL,)>*l PCS, S I GVAL.S I GPOS) 

IF  (NSTROKE(N) . EO. 1 )  THEN 
PYfMNVAL(J)  -  MWAL 
PWSD(J)  -  SIGVAL 
ELSE 

GEMNVAL(J)  -  MNVAL 
GESD(J)  -  SIGVAL 
ENO  IF 

END DO  I  ENO  OF  J  LOOP  (DIFFERENT  REGIONS) 

ENODO  I  END  OF  N  LCOP  (DIFFERENT  REVOLUTION  TYPES) 
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C 

C  •••  Writs  output  to  o  fils  with  format  massbol./  fils 
C 

ENCOOE  (40.’("[8LIVEN.DATA.MASS]MASS3AL.".I3)'.OUTFILE)  NRUN 
OPEN  (UNIT-10, FILE-OUTFILE.TYPE-'NEW’ .ACCESS-'SEQUENTIAL' . 

«  FORI*-  ’  FORMATTEn '  ) 
c 

writs  (10. ’(lx"  ")’) 

wri  te  (18. 1  (28x"  1  . . — —  |  ")') 

writs  (I0,'(28x"|  MASS8AL  RESULTS  |  ")’) 

writs  (18.'(28x"i  . . — |  ")’) 

writs  (10. '(U"  ")') 

writs  (10.’(28x"Oato  Sst:  ",o15)‘)infilsv(i) 
writs  (10, ’(lx"  ")') 
writs  (10.  ’  (lx”  ’  ’)’) 
c 

WRITE( 1 0 . ’ (5X ’ ’  POWER  STROKES  GAS 

S  EXCHANGE") ’) 

wri  ts(10,  ’  (5x"  -  - 

$ - ")’) 

WRITE(  1 0 .  ’  (5X '  ’REG  FRM  TO  MEAN  STD  DEV  MEAN 

$  STD  DEV")’) 

WR1TE( 1 8 . ’ (5X ’  ‘  f  (deg)  (deg)  (cu  mm)  (cu  mm)  (cu  nm) 

$  (cu  mm)  "  ) ’ ) 

wri  ts(l0,  ’(5x‘  ’ -  -  -  -  -  - 

$  - ")■) 

DO  L  -1 .  NLM3EG 

writs  (10. ’(lx"  ")’) 

WRITE  (10, ’ (5X. 12. 1X.F7.2.2X.F7.2.4(3X.F7.2)) ’ )L.  LWR8(L), 

$  UPR8(L) .  PWS^VAL(L).  PWSO(L).  GE)^VAL(L).  GESD(L) 

ENCDO 

writs  (10.  ’  (lx"  ’  ’ )  ' ) 
writs  (10,  ’  (lx"  ’  ’)  ’) 
writs  (  10  .  ’  (  1  x  "  ")’) 

writs  (10, ’(lx"  ")’) 

writs  ( 10. ‘ ( 1 x  "NOTE:  If  negative  crank  angls  regions  specified. 

$  thsn  POSITIVE  msan  values  indi-ate  ")’) 

writs  (10. ’(lx"  that  UPstroke  volumes  ars  GREATER  than 

$  DCWNstroks  volumes.  ") ’) 

writs  (10, ’(lx"  Standard  deviations  ars  always  pos i t i vs .  "  ) ’ ) 

ENOOO  I  END  OF  I  LOOP  (MULTIPLE  FILE) 

END  IF  I  END  OF  NCHOICE  5  IF  LOOP 
C 

C  •••  Final  Branch 
C 

5  CONTINUE 

IF  (NCHOICE. NE.0)  GOTO  10 
C 

999  CONTINUE 
C 

C 

C  Formats: 

C 

1000  FORMAT  (/’  ...  OIL  FILM  STATISTICAL  ANALYSIS  •••’.///. 

$T10,’(1)  Enter  files  to  analyze.  ’,/, 

$T10,’(2)  Sst  piston  regions.  ’,/, 

$T10,’(3)  Linear  Crankongls  Int erpo I  a t i on : ’ , 1 x , A3 ,/, 

$T 10, ’ ( 4)  Check  information  for  1  4  2.  ’ ,/, 

$T10,'(5)  Compute  ovsroges  and  write  to  massbcl  file.’,/, 

ST10, ' (0)  or  CTRLtZ  for  exit’///, 

$’  Enter  choice:  ’ ,$) 

1001  FORMAT  (/’  ...  OIL  FILM  STATISTICAL  ANALYSIS  »»*’,///. 

$T10.’(1)  Enter  files  to  analyze.  ’,/, 
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$T10,’(2)  Set  piston  regions.  ' ,/. 

$T10,'(3)  Linear  Crankongle  In t erpo I o t i on : ’ , 1 x , A3 ,/, 
$T10,‘(4)  Check  information  for  1  t  2.  * ./, 

$T10,‘(5)  Processing  file:  ’,A30,’  Be  patient!',/, 
CT”.  3 ,  ’  (0)  or  CTRLtZ  for  exit*///, 

$'  Enter  choice:  * ,$) 

STOP 

C 

END 
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F  .  3  "INTEGRATED:'.  FOR"  LIFTING 

C 

SUBROUTINE  INTEGRATEDF  (FI LM. LSD .UBD .VALUE. NRPTS . CRDL. CKRAD) 

C 

C  •••  This  suboutine  integrates  the  value  in  Film(2.i)  using  the  Crank- 
C  Nicholsen  method.  It  reads  the  step  size  from  the  first  row  of 

C  Film.  The  range  is  use r— spec i f i ed  with  lower  and  upper  bounds. 

0  The  difference  between  areas  of  successive  strokes  is  found  and 

C  returned  to  the  driver  program  for  averaging. 

C 

C  •••  DIMENSION  AND  DECLARE  VARIABLES 
C 

PARAMETER  NPMAX  »  250000  1  max  /  of  points/ch 

PARAMETER  NPCMAX  »  60000  !  max  /  of  points/cycle 

PARAMETER  NCMAX  -  250  I  max/  of  cycles 

C 

REAL  FI LM(2. NPCMAX) , LSD.UED.  VALUE(2) ,  AREA.  PCSA .  PCSB  . 

$  CRDL .  CKRAD .  DMED 
INTEGER  I.J.K.  NRPTS 
C 

C  •••  SET  INITIAL  CONDITIONS  BEFORE  INTEGRATING 
C 

J  -  1 
6  K  -  1 

DO  1-1. NRPTS 

IF  ( LEO . GE . FI LU( 1,1))  THEN 
K-I 

cc  a 

ccc  write(5,'(2f10.3)’)lbd,filn(l,k) 

ccc 

ELSE 
GOTO  3 
END  IF 
ENCCO 

S  CONTINUE 

c 

C  •  ••  Perform  integration 

V* 

AREA  -  0. 

10  IF  (FI LM( 1 ,K+1 ) . LE. USD)  THEN 

PCSA  -  CKRAD*CCS(FILM(1 ,K)/360.*2»3. H1593)+ 

JSC-RT (C.ROL*  *2  -  (CKRAD*  •  2 )  .  (SIN(FI  LM(  1  ,K)/360 .  *2*3 . 1 41 593) )  . *2) 

PCS3  -  CKRAD*CCS(FILM(1 ,K+1 )/3S3. *2*3. 14159)+ 

SSCRT (CRDL* • 2  -  (CKRAD* • 2) • (SIN(FI LV( 1 . X+1 )/350 . *2*3 . 1 41 59) ) • *2 ) 

C 

C  •••  Divide  by  1000.  to  convert  microns  to  mm 
C 

AREA-AREA+((FILM(2,K)+FILM(2,K-r1  ))/10B0.  )• 

$  A3S ( ( P0S3— PCSA)/2 . ) 

K  -  K+1 

ccc 

ccc  writa(6.'(3e10.3)’)poso,posb,orea 

ccc 

GOTO  10 
END  IF 

ccc 

ccc  write(5,’(2f10.3.e18.3)’)filra(1,k),film(2.k),area 

ccc 

C 

C  •••  Mirror  bounds  and  find  the  area  of  the  adjacent  downstroke, 

C  then  subtract  the  values  ond  return. 

C 

IF  (J.EQ.1)  THEN 
DMED  -  LSD 


O  O  O  ODD 
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LSD  -  -ueo 
USD  -  -OWED 
VAtUE(l)  -  AREA 
J  -  2 
GOTO  6 
ENO  [F 

VALU£( 1 )  -  VALUE(I)  -  AREA 
VALUE(2)  -  VALUE(I) 

Return  bounds  to  original  values 

DMBO  -  L90 
LBO  -  -UBD 
UBD  »  -CMED 


•  • •  Tha t ' s  the  end 

RETURN 

END 


o  o  o  o  o  o 
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F . 4  "AVERAGE. FOR"  LISTING 


SUBROUTINE  AVERAG E ( AVALS , MNVA L . MNPOS . S I GVA L.SIGPOS) 

...  SUBROUTINE  COMPUTES  AVERAGE  VALUES  AND  STANDARD  DEVIATION 
FOR  THE  FEATURE  OF  INTEREST  N  PASSED  IN  MATRIX  ’AVALS* 

...  DIMENSION  AND  DECLARE  VARIABLES 

PARAMETER  NPMAX  -  250000  I  max  f  of  points/ch 

PARAMETER  NPCMAX  -  60000  I  max  f  of  points/cycle 

PARAMETER  NCUAX  -  250  !  max#  of  cycles 

REAL  AVALS(2 .NCMAX) , MNVAL .MNPOS , SIGVAL. SIGPOS 
$  XSUM.  SUMPOS,  XVAL.  DUCT,  XDEV.  POSDEV,  SDSUM.  SDPSUM 
INTEGER  I.J.K 
1-1 

XSUL+0. 

sumpos-0 . 

10  IF  (AVALS (2 , I ) . NE. 0 . 0)  THEN 
XVAL  -  AVALS(2. I) 

XSUM  -  XSUM  +  XVAL 

SUMPOS  -  SUMPOS  +  AVALS(I.I) 

CUCT  -  I 
I  — 1  +  1 
GOTO  10 
END  IF 

UNVAL  -  XSLM/CMCT 

mnpcs  -  sumpcs/cm:t 

1-1 

SDSU+0. 

SDPSUV+0. 

12  I F ( AVALS ( 2 , I ) . NE . 0 ■ 0)  THEN 

XOEV  -  (AVALS(2 .  I  )  -*#FVAL)*»2 
POSDEV  -  (AVALS(I.I)  ~  fcHPCS)*»2 
SOSUM  -  SOSLAf  +  XOEV 
SDPSUM  -  SCPSLM  +  POSDEV 
CMCT  -  I 
I  — 1  +  1 
GOTO  12 
END  IF 

S I GVA  L-SCRT ( SDSUM/DUCT ) 

S I GPOS-SCRT ( SDPSUM/CMCT) 

RETURN 

END 
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Figure  1.  Engine  Instrumentation/Measurement  Equipment 
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Figure  2,  Piston  and  Rings  of  the  Kubota  EA300N 
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START 


Add  Radioactive 
Oil 


Operate  engine  for 
5  hours  to  burn  off 
light  ends 


Obtain  samples  in 
accordance  with 
Appendix  C. 


Prepare  samples  in 
accordance  with 
Appendix  D. 


Drain  Oil 
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oil 


I 
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flush  with 
flushing 
\  oil?  / 
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Figure  4.  Procedural  Flow  Chart 
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OIL  CONSUMPTION  USING  SAE-30 
(full  load  @  1500  and  3000  rpm) 
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RPM 


Figure  5 
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OIL  CONSUMPTION  USING  15W-40 
(full  load  @  1500  and  3000  rpm) 
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Figure  6 


GRAPH  OVERLAY  OF  SAE-30  AND  15W-40 
OIL  CONSUMPTION 
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Figure  7 


oil  consumption  (g/hr) 
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COMPARISON  OF  THE  OIL  CONSUMPTIONS 


Figure  8 
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VOLATILITY  OF  NEW  SAE-30 
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Figure  9 


(full  load  @  1500  rpm) 
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Figure  14.  Areas/Volumes  of  possible  oil 
consumption  mechanisms 
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LABELED  REGIONS  OF  COMP  STROKE: 
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Figure  15 
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LABELED  REGIONS  OF  COMP  STROKE 
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Figure  16 


OIL  FLOW  FOR  THE  POWER  EXCHANGE  STROKES 
USING  SAE-30  AT  1500  AND  3000  RPM 
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Figure  17 


OIL  FLOW  UP  THE  LINER  (g/hr) 
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OIL  FLOW  FOR  THE  POWER  EXCHANGE  STROKES 
USING  15W-40  AT  1500  AND  3000  RPM 
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Figure  18 


OIL  FLOW  UP  THE  LINER  (g/hr) 
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OIL  FLOW  FROM  THE  GAS  EXCHANGE  STROKES 
USING  15W-40  AT  1500  AND  3000  RPM 


PISTON  REGION 
Figure  20 
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OIL  FLOW  FOR  ALL  STROKES 
USING  SAE-30  AT  1500  AND  3000  RPM 


PISTON  REGION 
Figure  21 


OIL  FLOW  UP  THE  LINER  (g/hr) 
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OIL  FLOW  FOR  ALL  STROKES 
USING  15W-4C  AT  1500  AND  3000  RPM 


PISTON  REGION 


Figure  22 
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NON-DIMENSION ALIZED  OIL  CONSUMPTION 
AND  REGION  II  VALUES  FOR  POWER 
EXCHANGE  STROKES  USING  SAE-30 


Figure  23 


Ro'Ro'  (Ro’=  respective  value  at  1500  rpm) 
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NON-DIMENSIONALIZED  OIL  CONSUMPTION 
AND  REGION  II  VALUES  FOR  POWER 
EXCHANGE  STROKES  USING  15W-40 
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Figure  24 
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OIL  FLOW  RATES  OF  THE  UPPER 
HALF  OF  THE  PISTON 
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Figure  25 
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VISCOSITIES  OF  SAE-30  AND  15W-40 
vs.  RPM 


1000  2000  q  pm  3000  4000 


Figure  27 


surface  tension 
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SURFACE  TENSION  vs.RPM 


Figure  28 
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VISCOSITIES  vs.  OIL  CONSUMPTIONS 


Figure  29 


